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Abstract 


We demonstrate that the restoration of chiral symmetry at finite-T in a class 
of confining Dyson-Schwinger equation (DSE) models of QCD is a mean held 
transition, and that an accurate determination of the critical exponents using 
the chiral and thermal susceptibilities requires very small values of the current- 
quark mass: logio(m/mu) < —5. Other classes of DSE models characterised 
by qualitatively different interactions also exhibit a mean held transition. 
Incipient in this observation is the suggestion that mean held exponents are a 
result of the gap equation’s fermion substructure and not of the interaction. 
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It is anticipated that the restoration of chiral symmetry, which accompanies the formation 
of a qnark-glnon plasma at hnite temperatnre, T, is an eqnilibrinm, second-order phase 
transition. Snch transitions are completely characterised by two critical exponents: (/3,S), 
which describe the response of any one of the eqnivalent chiral order parameters, X, to 
changes in T and in the cnrrent-qnark mass, m. Denoting the critical temperature by Tc, 
and introducing the reduced-temperature t := T/Tc — 1 and reduced mass h := m/T, then 

X oc (—t —>• 0“ , h = 0 , (1) 

A' oc , h ^ 0+ , t = 0 . (2) 


Calculating the critical exponents is an important contemporary goal because of the widely 
conjectured notion of universality, which states that their values depend only on the symme¬ 
tries of the theory, the dimension of space, and whether or not the interaction is short-range. 
In this case many theories could be grouped into a single universality class labelled by (/d, 6), 
and chiral symmetry restoration in QCD might be describable by a much simpler model. 

The success of the nonlinear cr-model in describing long-wavelength pion dynamics un¬ 
derlies a conjecture |jl|] that the restoration of chiral symmetry at finite-T in 2-flavour QCD 
is characterised by a free energy whose interaction is just that of the 3-dimensional, = 4 
Heisenberg magnet (0(4) model). This theory has been studied extensively and its critical 
exponents are (3^ = 0.38 ± 0.01, 6^ = 4.82 ± 0.05. 

The conjecture that 2-flavour QCD is in the 0(4) universality class was reconsidered in 
Ref. 0, where it is argued that the compositeness of QCD’s mesons affects the nature of 
the phase transition. The Gross-Neveu model, dehned by the Euclidean action 


d 


d+1 


X 


V’(*7 -8 + 171 + g<j)il) — 


( 3 ) 


is presented as a counter-example. Arguments kindred to those of Ref. |I| suggest that the 
chiral transition in this model should be described by a d-dimensional Ising model because 
of the discrete a —> —a symmetry, which is manifest after integrating out the fermions. 

At T = 0 the Gross-Neveu model exhibits chiral symmetry restoration as the coupling 
g ^ Qc and at leading order in a 1/AA-expansion, where Np is the number of fermions, 
the critical exponents are easily calculated by exploiting the divergence that appears as 
the gap vanishes in the gap equation: /d-r=o = l/('^ ~ I); <5^=0 = T^e same analysis 
applies at £nite-T but in this case the absence of fermion Matsubara-frequency zero-modes 
ensures that the integrand in the gap equation is hnite as the gap vanishes, so that: ~ 

^MF _ [\/2), = (5^^ = 3. These exponents are apparently unchanged by l/Np- 

corrections in which case the hnite temperature transition is characterised by mean held 
exponents, independent of the dimension, not Ising model exponents. If this is correct then a 
description of the transition in terms of elementary meson helds is inadequate in the vicinity 
of the second-order phase transition 0. 

However, in a subsequent study |0 of a Yukawa model, which is equivalent to the Gross- 
Neveu model in the large-Yi? limit, it was argued that the chiral transition is described by 
the expected Ising model exponents but only on a domain of width ~ 3/Np in the vicinity 
of the critical-temperature and/or -external-held, with mean held exponents describing the 
evolution of the order parameters outside that domain. In this case 1/AA-correct ions do 
appear to be important. 
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These simple examples illustrate why the universality class containing QCD remains 
uncertain and why it will be difficult to determine if, as is plausible, the QCD chiral transition 
exhibits similar complexities; e.g., mean held exponents which may, or may not, evolve into 
0(4) exponents very near the critical-coupling and/or chiral limit. We explore this question 
by searching for a systematic trend in a comparison of {f3, 5) calculated in a range of classes 
of Dyson-Schwinger equation models of QCD. 

Calculating the critical exponents directly from Eqs. (p and (|^) is often difficult because 
of numerical noise near the critical temperature. Another method is to consider the chiral 
and thermal susceptibilities: 


Xh{t,h) 


d X(t, h) 
dh 


Xt{t, h) 


dX{t, h) 



(4) 


At each h, Xh{t, h) and Xt{t, h) are smooth functions of t with maxima at the pseudocritical 
points tpj. and where 

tjc oc oc . (5) 

Since j36 > 0, the pseudocritical points approach the critical point, t = 0, as h —0+ and 


Xr := Xhitpc, h) (xh Zh:=l-^, (6) 

xr ■= Xt{tl^, h) oc , zt:= p). (7) 

Therefore, by locating the pseudocritical points and plotting the peak-height of the suscep¬ 
tibilities as a function of h one can determine Tc, (3 and 5. 

This method was employed in numerical simulations of finite-T lattice-QCD on an 
8^ X 4-lattice with three values of the current-quark mass: = 0.075,0.0375,0.02, 

which for 1/a ~ 4Tc, Tc ~ 150 MeV, correspond to = 45,23,12 MeV, and yielded: 
plat _ g gg g gg^ ^lat _ ^ g g rjpg values are consistent with a straightforward 
realisation of the 0(4) hypothesis, however, the errors are large. The situation deteriorates 
further when more recent simulations on larger lattices and with smaller quark masses are 
taken into account 0. A value of 1/h*®'* = 0 is consistent with the new data, and this can 
be interpreted as a signal of a first order transition. Our calculations indicate that this 
unexpected result may be an artefact, due in part to the large values of the quark mass that 
contemporary lattice-QCD simulations are restricted to. 

Herein we analyse Xh{t,h) and Xt{t,h) in a class of conhning DSE models that under¬ 
lies many successful phenomenological applications at both zero and £nite-(T,/i) [§]. The 
foundation of our study, as for many others PH43|, is the renormalised quark DSE 

■= ■ P A{P^k) + ix^u^kCip^^) + B{p^^) (8) 

= Z2ix-P + Z2 {iXAi^k + rn^m) + S'(p^J. (9) 

Here (Pi^*,) := {p,uJk) with ujk = (2fc -|- 1) vrT the fermion Matsubara frequency, and mbm is 
the Lagrangian current-quark bare mass. The regularised self energy is 


^'ip<^k) = W ■ P^AiPco^) + ^74 Wfc S^(p^J + , 

/*A 

^'APu^k) = / tAD^^{p-q,uJk-uJi)ltT[VrX^S{qu;,)T^{q^PP^^)] , 

Jl,q 
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( 10 ) 

( 11 ) 






where: T = A,B,C] A,B,C are functions only of and a;|; Va ■= 

Vb '■= ■Z’l, Vc '■= —(2'i/a;fc)i74; and := T J2^-oo with a mnemonic to 

represent a translationally invariant regularisation of the integral and A the regularisation 
mass-scale. In renormalising the quark DSE we require 




\p\^+u;^0=e 


ii-p + *74 ujo + mn , 


( 12 ) 


and hence the renormalised self energies are 


r{p.,;0 = ^^ + si■(p...;A)-Ei-(C-;A). 


(13) 


(Co)^ := ^ ^ AB,C, (a ^ 1 =&, and {b = mK((). 


Equations define the exact QCD ffap equation. It has obvious qualitative simi- 

larities to the Gross-Neveu gap equation, in particular the explicit dependence on ujk- 

In Eq. (pd]) , is the renormalised dressed-quark-gluon vertex whose structure 

at T = 0 has been much considered 0- It is a connected, irreducible three-point function 
that should not exhibit light-cone singularities in covariant gauges; i.e., it should be regular 
at {p — ^‘^ + {uJk — uJi)'^ = 0 1^. A number of Ansdtze with this property have been proposed 
and employed, and it has become clear that the judicious use of the rainbow truncation 
= Ju) in Landau gauge provides phenomenologically reliable results [|I^. We 
employ it herein, in which case a mutually consistent constraint is Zi = and Z^ = Z^- 
The rainbow truncation is the leading term in a 1/iVc-expansion of Ty{q^^pp^^). 

D/MuiPn^) in Eq. (Plj) is the renormalised dressed-gluon propagator, which has the form 






fiiy \ 


(14) 


:= 


0, p and/or z/ = 4, 

r PiPj • ■ 1 O Q 

dij - p,u = t,j = 1,2,3 , 

pZ 


(15) 


with PjLiPnJ + P^^ipnj = ion - PuPulT,i.i PaPa\ P,o = 1,...,4. (Qt := 2wkT is the 
boson Matsubara frequency.) The primary class of DSE models we consider is that in which 


We write 


the long-range part of the interaction, g^Dniu{pQ^,), is an integrable infrared singularity 


Af(poJ = V{pnp,mg), l^oipn^) =T>(po,;0), 
V{pQp,mg) := 2TT^D^5okS^{i^ + T>m(po;,; m^) , 


(16) 

(17) 


where D is a mass-scale parameter and T>M{pnkimg) may be large in the vicinity of p^^ = 0 
but must be hnite. This type of infrared singularity is motivated by T = 0 studies of the 
gluon DSE. Studies in axial gauge |0, where ghost contributions are absent, and in Landau 
gauge [0, when their contributions are small, indicate that Dgy{k) is signihcantly enhanced 


in the vicinity of /c^ = 0 relative to a free gauge-boson propagator. In the neighbourhood 
of = 0 the solution is a regularisation of l//c^ as a distribution, and the enhancement 
persists to ~ 1 GeV^ where a perturbative analysis becomes quantitatively reliable. 

The model of Refs. is obtained with D := 7^/2, 'DM-.=A{pnA'^g) = 0’ mass- 


scale p = 1.06 GeV hxed [0 by htting tt- and p-meson masses at T = 0. It is an ultraviolet 
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finite model and hence the renormalisation point and cutoff can be removed simultaneously. 
A current-quark mass of m = 12 MeV yields tUtt = 140 MeV. As an infrared-dominant model, 
it represents the behaviour of poorly away fromp^^ 0. However, the artefacts 

this introduces are easily identified and the model exhibits many features in common with 
more sophisticated Ansdtze. 

is obtained with D := (8/9) and 


The model of Ref. 


'I^M:=B{Pnk]'rng) = 


1 _ 




(18) 


SUk •= Poj, where = (8/3) vr^T^. The mass-scale = 0.69 GeV = 1/0.29 fm marks 
the boundary between the perturbative and nonperturbative domains, and was also hxed 


by requiring a good description of vr- and p-meson properties at T = 0. At a renormalisation 
point of C = 9.47 GeV, niR = 1.1 MeV yields 771^, = 140 MeV. This model adds a Goulomb- 
like short-range interaction to that of Refs. 0 , thereby improving its ultraviolet behaviour. 
The ratio of the coefficients in the two terms of Eq. (O) is chosen so that the long-range 
effects associated with SokS^{p) are completely cancelled at short-distances [p0| . 

We also consider a £nite-T extension of Ref. [M dehned by 




'DM-.=c{Pn,;rng) =—Dsn,e 




In 


r (l-f snk/A-QCD 


g-sr!fc/(4m?) 

SOi. 


(19) 


with r = — 1, 7m = 12/25, = (16/5) tt^T^, and = 234 MeV. This further 

improves the ultraviolet behaviour, via the inclusion of the one-loop In-suppression at sqj, S> 
Aqcd, and through the Gaussian incorporates the intermediate-range enhancement observed 
in zero-temperature gluon DSE studies. The parameters D, rrit and oj can be fixed at 
T = 0 by requiring a good £t to a range of vr- and A'-meson properties. Herein we consider 
two phenomenologically equivalent parameter sets that differ only in the value of a;: D = 
0.78 GeV^, rrit = 0.5 GeV, and: lUJ = 0.6 mj, 2 ^ = 1.2 mt. At T = 0 in this model the quark 
mass function, M{p^) := B[p"^)/A{p‘^), evolves according to the one-loop renormalisation 


group formula for p^ > 


20 Aq(-;j3, 


and with lOO a renormalisation point invariant current- 


quark mass of m = 6.6 MeV yields = 140 MeV while with 2 <^, rh = 5.7 MeV effects 
that. 

An often used order parameter for dynamical chiral symmetry breaking is the quark 
condensate 


16 


- (??)c := Nc Jim ^4(C, A) [ tr^ [5'o(Pa;,)] , (20) 

A^OO Jl,p 

for each massless quark flavour. Here the subscript “0” denotes that the dressed-quark 
propagator is a chiral limit solution of Eq. (^, and Z 4 ((',A) is the mass renormalisation 
constant: Z 4 ((/, A) mij((/) = Z 2 (C, A) mbm(A). There are other, equivalent order parameters 
and from Eq. (|20|) it follows that one such is X := Bo{p = 0, cuq). This being a bona fide order 
parameter is particularly useful and important because it means that the lowest Matsubara 
frequency completely determines the character of the chiral phase transition. 
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The quark DSE obtained with Va is an algebraic equation. Its solution is therefore 
easy to analyse and either directly, via Eqs. o and (H), or using the susceptibilities, it 
is straightforward to establish P that this model has mean held critical exponents and to 
determine the critical temperature in Table |. The exponents are unchanged (e.g., ^"(0, is 
linear on h G 10“®]) and Tc reduced by < 2% (Tc = rj/{2 ti) = 0.15915 r/ —> 0.15646 r/) 

upon the inclusion of the higher-order l/A');-corrections to the dressed-quark-gluon vertex 
discussed explicitly in Refs. [PH] . 

The critical behaviour of the "D^-model is harder to explore because the solution of the 
quark DSE must be obtained numerically. The apparently simple interaction is actually more 
complicated to study than that dehned with Vc because the ultraviolet In-suppression 
characteristic of asymptotically free theories is lacking. Therefore the renormalisation group 
properties are those of quenched QED, which has notable pathologies . 

Chiral symmetry restoration in this model was studied in Ref. ^j. Herein, however, 
we report results obtained with an improved numerical procedure. The most signihcant 
change, implemented in order to make recovery of the T = 0 limit easier, was to place a 
hxed ultraviolet cutoff, 4 A, on := q'^+uf. This requires that a different three-momentum 
cutoff and grid be used for each Matsubara frequency. For small-T this corresponds closely 
to an 0(4) invariant cutoff procedure, and for all T it provides for an accurate discrete 
representation of the kernel. In the calculations we chose 4 A ~ 1-5 C ~ 14GeV, and for 
a given value of T the number of Matsubara modes required, /max, is determined from 


~ 4A. 

We employ two chiral order parameters in analysing chiral symmetry restoration 


X B{p — 0,a;o), Xq : — 


B{p = 0,a;o) 
0(p = 0,a;o)' 


( 21 ) 


They should be equivalent and, as we will see, the onset of that equivalence is a good 
way to determine the h-domain on which Eqs. ®-( 0 ) are valid. Further, we have verihed 
numerically that for m = 0 and f ~ 0 : /jr oc {qq) oc ^(^, 0 ); i.e., that these quantities are 
all bona fide order parameters. It thus follows from the pseudoscalar mass formula |]T 6 
f^ml = 2 m/j(C)(gg)°, that m. 


diverges at the critical temperature 0 , 


m„ 


oc i-t) 


-0 


t O' 


( 22 ) 


The critical temperature in Table H was obtained^ using Eqs. (^). Its value is insensitive to 
whether or is used and also to the order parameter. Xh^W and xf‘^(h), obtained using 
both order parameters, are depicted in Fig. ^ and, following Eqs. (H) and (|^, the critical 
exponents are determined from these curves. In comparison with Ref. [Q it is obvious that 
there is signihcantly less error in the function values. This makes it possible to distinguish 
the nonzero curvature on the domain —4.5 < logig h < —3, whereas in Ref. ^ the error in 
the function values was such that the scaling relations, Eqs. (i-®. appeared to be valid on 
this domain. In fact, the scaling relations are not valid until 


^ In these calculations we corrected a numerical error in Refs. PJT^: the Debye mass was a factor 
of three too large therein. That is the origin of the 15% increase in Tc, which also serves to illustrate 
the small influence of the perturbative Debye mass on the results. 
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Iogio(ft/A«) < -7. 


( 23 ) 


hu := mji/Tc is defined with the current-quark mass that gives = 140 MeV in this model. 

The scaling domain is most easily identihed by dehning a “local” critical exponent for 
each of the equivalent chiral order parameters: 


inxr-inxr:i 

In hi — In hi+i ’ 


(24) 


where and are adjacent data pairs. Zi{h) for the present model is 

depicted in Fig. |^. h lies in the scaling region when Zi is independent of the order parameter. 
Anticipating the curvature, we £t the susceptibilities using the functional form 


logio X’’" = 


1 -h ao ? ai 


(26) 


? = logio ^1 over the smallest hve values of and the critical exponent is then: 2 ; = 01 / 03 . 
We estimate the error by systematically eliminating one value of ? at a time from the fit 
and observing the change in z. In this way we obtain the values of Zh and Zt in Table |I[ 

A linear fit on G [—4.5, —3] yields: Zh = 0.79, Zt = 0.40, which may be compared 
with the values in Ref. P: Zh = 0.77 ± 0.02, Zt = 0.28 ± 0.04. Comparing the values of Zh 
highlights the error introduced by incorrectly identifying the scaling region. The comparison 
between the values of Zt indicates our significant improvement of the numerical method. 

The study of chiral symmetry restoration in the model obtained with Vc is straight¬ 
forward, with the additional In[sn*]-suppression in the ultraviolet making the numerical 
analysis much simpler. The critical temperature for each parameter set is presented in Ta¬ 
ble in, and on this cn-domain: Tdui) ~ 88 -|- 107 ui. xTW Xt^{h) for the smallest value of 
CO, which corresponds to the parameter set in Ref. [|T^], are depicted in Fig. |^. The results 
for the larger value of a; are similar. In these cases the scaling relations are only valid for 


icn: logio(h/h„) < -5, 20 ;: logio(h/h„) < - 6 . (26) 

It is clear from Table | that each of these models is mean held in nature. In hindsight 
that may appear unsurprising because the long-range part of the interaction is identical 
in each case and the correlation length diverges as t —>■ 0. However, the models differ in 
detail by the manner in which the interactions approach their long-range limits, and our 
numerical demonstration of their equivalence required a very careful analysis and extremely 
small values of the current-quark mass, Eqs. (|23|) and (p^). 

This last observation is also likely to be true in QCD; i.e., while the critical temperature is 
relatively easy to determine, very small current-quark masses may be necessary to accurately 
calculate the critical exponents from the chiral and thermal susceptibilities. If that is the 
case, the calculation of these exponents via numerical simulations of lattice-QCD will not 
be feasible. The discrepancies described in Ref. 0 could be a signal of this. 

The class of models we have considered explicitly preserves the chiral Ward-Takahashi 
identity and can describe the long-wavelength dynamics of QCD very well |]8|JT^. Impor¬ 
tantly, it describes that dynamics in terms of mesons that are quark-antiquark composites. 
The characteristic feature of this class is the behaviour of the long-range, conhning interac¬ 
tion. It provides an additive, algebraic driving term in the quark DSE that is proportional 
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to the dressed-quark propagator, which means that boson Matsubara-frequency zero-modes 
do not influence the critical behaviour determined from the gap equation. 

The class of Coulomb gauge models considered in Ref. 0 also describes mesons as com¬ 
posite particles and it too exhibits mean held critical exponents. The long-range part of the 
interaction in that class of models corresponds directly to the regularised Fourier amplitude 
of a linearly rising potential. Hence it is not equivalent to ours in any simple way, except 
insofar as zero modes do not influence the gap equation. 

0 


We have also considered the class of models obtained with 


'^{Pujk - Qujpmg) oc 
This class includes 


g{\^), where g{\p\) is a non-increasing function of its argument, 
the non-confining NJL-model, which is obtained with g{\f^) = 0{1 — Ipj/aA), where 3 A. is a 
cutoff parameter, and can also provide a good description of long-wavelength pion dynamics. 
The models in this class describe mesons as composites and exhibit an explicit fermion 
substructure, and those we have explored have mean held critical exponents. The same is 
true of separable models with gi = gtiuil +p^) P3| . 

The quark DSE is the QCD gap equation and the many equivalent chiral order parameters 
are directly related to properties of its solution. We have observed that four classes of 
models exhibit the same (mean held) critical exponents. The classes are distinguished by 
the qualitatively distinct type of interaction they employ in the gap equation but similar in 
using the rainbow truncation. Only in our simplest confining model did we consider the effect 
of 1 /W-corrections to the quark-gluon vertex, and in that case the critical exponents were 
unchanged. Incipient in these results is the suggestion that mean held exponents are a feature 
of the essential fermion substructure in the gap equation, and not of the interaction.^ If that 
is correct then chiral symmetry restoration at flnite-T in QCD is a mean held transition. 
This hypothesis is supported by chiral random matrix models of QCD, whose formulation 
relies only on the symmetries of the Dirac operator and which also yield mean held 
exponents. It can likely only be false if 1/iVc-corrections to the vertex are large in the 
vicinity of the transition, in which case their null effect in our simplest model will have been 
misleading. 
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TABLES 



A 

B 

Cl.; 


Tc{MeV) 

169 

174 

120 

152 

Zh 

0.666 

0.67 ± 0.01 

0.667 ± 0.001 

0.669 ± 0.005 

Zt 

0.335 

0.33 ± 0.02 

0.333 ± 0.001 

0.33 ± 0.01 


TABLE I. Critical temperature for chiral symmetry restoration and critical exponents charac¬ 
terising the second-order transition in the four exemplary models. Mean field critical exponents 
are: Zh = 2/3, zt = 1/3. 
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FIGURES 



h 

FIG. 1. y^{h) (circles) and (squares) calculated in the model defined via Vb in Eq. (p^: 

filled symbols - X, open symbols - Xc- The slope of the straight lines is given in Table | and they 
are drawn through the two smallest h-values. 
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FIG. 2. (circles) and z\ (squares) from Eq. ( p^ ) in the model defined via Vq in Eq. ( pSj ) 
filled symbols - open symbols - Xq- The dashed lines are the mean field values: Zh = 2/3 
= 1/3. 
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FIG. 3. X^^{h) (circles) and Xt^W (squares) calculated in the model defined via Vc in Eq. (W) 
with CO = 0.6 m*: filled symbols - X, open symbols - Xc- The slope of the straight lines is given in 
Table | and they are drawn through the two smallest h-values. 
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